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Minimal charge gap in the Ionic Hubbard Model 
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We study the ionic Hubbard model at temperature T = within the mean-field approximation and 
show that the charge gap does not close completely at the ionic-band insulator to antiferromagnetic 
insulator transition, contrary to previous expectations. Furthermore, we find a new intermediate 
phase for on-site repulsions U > Uc for different lattices and calculate the phase diagram for the 
ionic Hubbard model with alternating U, corresponding to a Cu-O lattice. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The generalization of the Hubbard model with different 
on-site energies on the neighboring sites has been named 
the ionic Hubbard model (IHM), 



H 



i.( 
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(1) 
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which is characterized by a nearest-neighbor hopping am- 
plitude t, an on-site Coulomb-repulsion U and a site- 
dependent on-site energy ±£'o. Here |i| =evcn/odd re- 
spectively for a A- and a B- site on a bipartite lattice. 

The IHM is used in two contexts: (a) For the descrip- 
tion of the neutral-ionic transition (NIT) in organic 
mixed-stack charge-transfer (CT) crystals. The stacks in 
CT crystals form quasi-ld insulating chains with alter- 
nating donor and acceptor molecules. The charge on the 
acceptor p and on the donor —p characterize the crystal 
state. For p < 0.5 the crystal is said to be neutral-like, 
otherwise it is ion-like. The transition from one region 
into the other due to the change of temperature or pres- 
sure is called NIT. (b) The IHM has also been used in 
the context of ferroelectrics and superconductivityin the 
transition-metal oxides 01 . It has been argued ]2| that 
the influence of the underlying lattice on the electronic 
system could be large in the critical region, leading to a 
non-linear electronic polarizability. 

The phase diagram of Q has been discussed previ- 
ously by several authors. Oritz et al. did find a sin- 
gle phase transition in a mean-field study using a one- 
dimensional density of states The transition takes 
place between paramagnetic and the antiferromagnetic 
state. The point where it takes place has a special role 
with respect to the electrical conductivity. Whereas the 
system is a semimetal in this point, in other points of 
the phase diagram it is an insulator. The possibility of 



a semimetallic transition point involving antiferromag- 
netism has also been discussed in a LDA-study H. 

Gidopoulos et al. performed a mean-field decoupling 
for two-dimensional bipartite square and honeycomb lat- 
tices 0. They did find two special values of U for a 
given energy-alternation Eq . The first denotes the phase- 
transition where the magnetic order sets in and at the sec- 
ond special value for the one-site repulsion U the charge- 
gap for one spin species would presumably close. In what 
follows, we'll show that there arc actually two phase tran- 
sitions in the mean-field decoupling scheme. Further- 
more, it will be shown that the solution where the gap 
for one spin species vanishes is thermodynamically not 
stable, the charge-gap consequently does not close at any 
point in the phase diagram. 
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FIG. f: Phase diagram for the Id DOS. For U < Uc the 
intermediate phase (I) is thermodynamically not stable and 
the transition from antiferromagnetic (AF-I-CO) to paramag- 
netic phase (CO) is of the 1st order. The U > Uc region has 
a stable (I)-phase. The transition between the two magnetic 
phases (AF-I-CO and I) is of the 1st order, the (I)-phase goes 
into the (CO)-phase through the 2nd order phase transition. 
The inset shows asymptotic behavior of the phase border for 
small U. 
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II. ORDER PARAMETERS 

Application of the mean-field decoupling 

leaves us with the system of equations 

oo 

Eo + UA-a 



(2) 



Aa = - deD{e) 



v/e2 + {Ea + C/A_,)2 



(3) 



for the two order parameters Ao- defined through the re- 
spective spin-densities rij^/^^ = 1/2 =p on A- and B- 
sites. D{e) is a free-particle density of states (DOS). In 
what follows, Id-DOS will be used, but qualitatively the 
same results have been obtained for the flat and the semi- 
circular DOS. Quasi-particle states are accommodated in 
four bands with dispersion 



U 



a = ±l. (4) 



In the ground-state the a = —1 bands are filled. The 
thermodynamic stability of a possible self-consistent so- 
lution of ^ is determined, at T = 0, by the total energy 
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Above equations are valid for the half-filled system. Use 
of the particle-hole symmetry 



c,.\FB) = (-1) 



|j+i|wt 



|0). 



(6) 



where \FB) denotes a full band and |0)h a hole vacuum 
0, shows that ep = U /2. Together with the dispersion 
of the bands, this implies that the number of each spin 
species is the same. 



III. PHASE DIAGRAM 

The phase diagram of the IHM is shown on Fig. ^ 
Regarding the [/-axis, it has two separate regions. For 
U < Uc ^ 4.25i one finds two phases: the charge-ordered 
(CO)-phase and a mixture of the charge ordered and the 
antiferromagnetic phase (AF-I-CO). In the U > Uc re- 
gion there are three phases. Beside the (CO) and the 
(AF-f CO)-phases, we find another phase which we will 
call intermediate (I)-phase. The (I)-phase has the same 
order parameters as (AF-l-CO)-phase. 

The phases shown in Fig. ^ can be observed clearly 
in the curves for the staggered magnetization nist = 
|A| — A|| shown in Fig. |21as a function of Eq and var- 
ious values of U. We notice a jump in the curves for 




FIG. 2: Staggered magnetization curves for diflterent [/'s. 
Blank, inclosed region contains jumps of the ruat- Shaded 
regions are placed under the magnetization curves of the (I)- 
phase. 



all values of U, indicative for the 1st order phase tran- 
sition. The jumps define a region of m^t for which the 
system is unstable, illustrated by the light shaded area 
in Fig. 121 When U > Uc there is no direct transition 
from the Iow-Eq (AF4-C0)-phase to the high-£;o (CO)- 
phase which is characterized by rUst — 0. The (AF-I-CO)- 
phase now shares a border with (I)-phase. The jump in 
the staggered magnetization vanishes asymptotically as 
U becomes large. The transition from (I) to (CO) is, on 
the other hand, of the 2nd order. From the magnetiza- 
tion curves it is also transparent that the width of (I) 
rises as one increases U (illustrated by the dark-shaded 
areas in Fig. |2J). 

The first-order nature of the transition between the 
(AF-l-CO)-phase realized for low values of Eq and the 
(CO) and (I)-phase respectively shows up in a jump of 
the charge-gap, as shown in Fig. O The discontinuity in 
the charge-gap is manifested through two different curves 
of the gap as a function of U for the two-sides of the 
transition respectively. 



(AF+CO) side 
(CO) and (I) side 





FIG. 3: Charge gap as a function of U on the 1st order tran- 
sition point (AF-(-CO)^(I) for U > Uc and (AF-HCO)^(CO) 
for U < Uc- On the left plot we zoomed in the curves show- 
ing behavior of the gaps in the region of the parameter space 
where the {AF + CO) phase shrinks rapidly. 
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A difFcrence between (I) and (AF+CO) is in the behav- 
ior of the charge gap as a function of E^, see Fig.0| The 
transition point between the (AF+CO) and (I)-phase is 
characterized by a minimum in the charge gap, see Fig.|21 
Even though the gap is close to zero, it remains finite. 
In the (AF+CO)-phase (Fig. the charge gap, as a 
function of Eq, decreases while it increases in the (I)- 
phasc (Fig. ^). As expected, this increase proceeds in 
the (CO)-phasc (Fig.H;). 
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FIG. 4: Evolution of the gaps for each spin direction as a 
function of Eq, for fixed U = 4\/2 t. The smaller gap is 
equivalent to the charge gap. (a) (AF+CO)-phase, the charge 
gap decreases, (b) (I) state, the charge gap increases, (c) 
(CO)-phase, the gap for both spin directions is the same. 



IV. THERMODYNAMIC STABILITY 

A self-consistent solution of Eq. (j2J) with finite values 
for the order parameters is not a guarantee that the 
corresponding phase is indeed realized in the system. The 
condition for the stability of a phase on T = is that 
it has the minimal energy with respect to the energies 
of other phases. This fact is important in the case of 
(AF-hCO)-(I) transition. There, the (AF-hCO) and the 
(I)-solutions show a hysteresis effect and the ground state 
is determined on the basis of the minimum energy prin- 
ciple. 

The examination of the energy space for the case 
U < Uc is shown on Fig.lSJa). On it one sees the energy 
Eq. (jSJ for the four different solutions which solve the self- 
consistency Eq. Q for the order parameters Ao-. They 
arc the three solutions (AF-|-CO,line), (CO, large filled 
circles) and intermediate phase (I, small filled circles) 
which we have already discussed, together with the fourth 
solution, which is always unstable (empty squares). The 
energy of the (I)-phase is larger than the one of the 
(AF-l-CO)-phase in the whole interval of its existence. 
(I)-phase is thus unstable. Crossing of the energies of 
(AF-I-CO) and (CO)-phase discloses an underlying 1st 
order phase transition. The lower inset shows energies of 
all phases in the whole Eq interval. 



For U = 3.2^/2 t > Uc the results are given in Fig.jSfb). 
Part of the (I)-phase is now stable. 

Stable solutions have been obtained by iterating 
Eq. directly and the unstable solutions by fixing one 
of Act's and allowing Eq to change. The fixed order- 
parameter Act has been chosen such that it is in the range 
which is not covered by the stable solutions. We have 
also tried to use fixed U , but the self-consistency map, 
for the parameters we examined, didn't turn out to be 
attractive and thus useful. Limiting cases: (1) Eq = Q and 
(2) J7 = assure us that in this limits we found all solu- 
tions. Thus, possible solutions with the weakly attractive 
self-consistency map are a possibility in the intermediate 
range of parameters. We cannot exclude possibility of 
their existence, but due to a few self-consistency arrange- 
ments we find them improbable. 

The charge gap vanishes when Act = —Eq/U for one 
spin species and A_ct = for the other species. The en- 
ergy of this solution is given in Fig. [SJa) and Fig. [Sfb) 
and is denoted by an empty triangle. It turned out to 
be unstable in the calculations done. Furthermore, from 
the performed calculations (see Fig. (S)) we see that the 
vanishing-gap solution has always (AF-I-CO) and (CO) 
companions with the smaller energy. This leads us to 
extend the conclusions obtained from our results to the 
whole parameter region, including interval U/t < 0.5 
which we haven't investigated numerically. Thus, the 
charge gap closes only in the point Eq ~ U = 0. 

Finally, wc would like to mention that the three sta- 
ble phases (CO), (I) and (CO+AF) arc characterized by 
different distribution of the charge densities Uj^/bo- = 
1/2 =F Act on A- and B-sites. We find: A| = A^ for 
the (CO) phase, A| ^ A^ and A|A| > for the (I) 
phase and A^Aj < for the (AF+CO) phase. 

V. ALTERNATING U IONIC HUBBARD 
MODEL (AIHM) 

A natural generalization of the ionic Hubbard model 
is the model where on-site Coulomb interaction on the 
atoms A and B is not the same. The Hamiltonian of the 
new system may be written as 

Haihm = HiHM + ^fi[/(-l)l*lni|nix (7) 

i 

For simplicity we consider here the extreme case where 
the Coulomb interaction disappears on site A. This im- 
plies dU ~ —U and a value of 2U for the Coulomb re- 
pulsion on site B. This model is than equivalent to the 
Cu-0 lattice model Q with correlated S-sitcs (Copper, 
lower on-site energy) and uncorrelated ^-sites (Oxygen, 
higher on-site energy). 

Let's write Eq = U/2 + cq and fi = U/2 + ^' , where /i 
is the chemical potential for the half-filled system. The 
transformation 

Eo = U/2 + eo ^ E'q^ U/2 - cq 
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FIG. 5: (a) Energy of the various states as a function of Eg. U = 2.9\/2 t < Uc is chosen such that the (I)-phase is unstable. 
The main panel is blowup of the energies shown in the inset, transformed by the subtraction of a suitably chosen straight so 
that the particular transitions can be easier visualized. The large empty triangle denotes the position of the vanishing charge 
gap in the thermodynamically unstable solution, (b) U = 3.2\/2 t > Uc so that the part of (I)-phase is stable. 



and 

fi = U/2 + fi' ^ = U/2 - ^l' 

leaves us with an equivalent Hamiltonian. This can be 
shown by performing the canonical particle-hole trans- 
formation Cj^cr (^1)''''^'] -cr- It yields 

<i,j>,a i 

"E(-1)''''(^+m'K. + 2A^m' (8) 

The difference with respect to the original Hamiltonian 
O lies in the constant term 2Nfi'. This, together with 
the fact that the order parameter operators preserve the 
same form in the hole picture, implies that the phase 
diagram of Q is symmetric with respect to the line Eq = 
U/2 in the parameter Eq and for dU = —U. 

The self-consistency equation, derived under the as- 
sumption that the ferromagnetic order parameter van- 
ishes, is given by 

oo 

where = Ea + C/A_„ + dU/2. 

The energy spectrum has a form A±.o-(e) = Aa- ± 
^/e^ + Bl where = U/2 + dUA-„. 

In Fig. we present the phase diagram of the AIHM. 
It contains the same phases as the IHM phase diagram. 
The dash-dot line in Fig. |3 indicates the symmetry axis. 
On it the system is antiferromagnetic with the vanish- 
ing charge-order parameter (AF). In the region around 



the symmetry line the system is in the (AF-I-CO) phase 
which, for U < Uc, makes transition into the pure charge- 
ordered state. For U > Uc the transition takes place into 
the (I)-phase and than into the (CO)-phase. 
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FIG. 6: Phase diagram of the half-filled AIHM with the Id 
density-of-states. It contains three phases. The charge or- 
dered (CO)-phase is given by a white area. The dotted area 
represents a mixture of the antiferromagnetic and the charged 
phase (AF-I-CO) and the black region covers a parameter 
range where intermediate (I)-phase appears. The symmetry 
of the phase diagram is denoted by the dash-dot line. 

Let's define a singlet gap as a minimal energy needed 
for a transfer of the electron in an empty state with- 
out a spin-flip and the triplet gap as a minimal energy 
of the transition where the spin is flipped and the total 
5"^ changes. In the mean-field formulation this gaps are 
given by A, = 2\B,\ and At = \A^ + \B^ \ + \B^ \ - 

A distinct property of the AIHM is the closing of the 
triplet gap in the (AF-I-CO) phase. This is illustrated 
in Fig. 13 IHM case with a uniform U is shown on 
Fig.EJa). In Fig.CJc) we see that the triplet gap vanishes 
for Ua — and Ub — 2U. A necessary (but not suffi- 
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FIG. 7: Singlet (dashed curve) and triplet (full curve) gaps 
as a function of the on-site chemical potential amplitude iJo 
for [/ = 6t. (a) JJ is the same for both atom species, (b) 
U A < Ub (c) U on site A vanishes. 

cient) condition for the vanishing of the triplet gap can be 
derived from Eq.0 and expressed in the form A|A| < 0. 
This is the case in the (CO+AF) phase. A|A^ > and 
the triplet gap is non- vanishing as can be seen in Fig. \7\ 
The exception to the condition A|A| < for disappear- 
ance of the triplet gap is a solution A| = A| = which 
exists for Eq = 

VI. DISCUSSION 

Previous mean-field studies of the ionic Hubbard model 
did find a vanishing charge gap at the transition point be- 
tween antiferromagnetic and paramagnetic state in 
contrast to our result of a minimal charge gap. Presented 
data has been obtained for a Id one-electron density-of- 



state. For comparison we have carried out calculation 
also for the flat and the semicircular density-of-states. 
The results changed only quantitatively, all major fea- 
tures discussed previously remain valid. 

It is interesting to compare with the Id-ionic Hubbard 
model, which was studied by Resta and Sorella Ul us- 
ing the boundary-condition- integration technique [g , by 
Brune et al. using DMRG and by Wilkins and Martin 
using QMC Fehske et al. considered the dynamical 
IHM coupled to phonons [llj . 

Fabrizio et al., in a Bosonization-study proposed 
the existence of a dimerized intermediate phase in the 
Id ionic Hubbard model. DMRG- and QMC-results 
indicate |2, that this bond-ordered (dimerized) 

insulating-phase extends to arbitrarily large values 
of U, due to the instability of the Id Mott-Hubbard 
insulator towards dimerization. DMRG studies of the 
Id-model 0,01, found indications of a 'strange' metalhc 
point with a finite charge gap at the transition from 
the band-insulator to the correlated insulating-state. 
We may speculate, that our mean-field result of two 
distinct phase transitions reflects, on the other hand, 
the occurance of magnetic long-range order possible in 
dimensions larger than one. We note, in this context, 
that the intermediate phase (I) evidenced in Fig. ^ and 
Fig. 1^1 does not show spontaneous dimerization. 



VII. CONCLUDING REMARKS 

The mean-field decoupling scheme is an approximation 
to the full evaluation of the Hamiltonian. Its validity 
is restricted to the small U limit and large dimension 
where it is qualitatively correct Q. We have found, 
however, that the stability-analysis for the various so- 
lutions possible for the ionic Hubbard model is highly 
non-trivial. Here we have presented two new results: (i) 
A non-vanishing charge gap for all parameters and (ii) 
the existence of two distinct phase transitions for larger 
values of U. 
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